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Abstract 

We have studied charged nuclear fragments produced by 200 - 400 MeV/nucleon carbon ions, 
interacting with water and polycarbonate, using a newly developed emulsion detector. Total and 
partial charge-changing cross sections for the production of B, Be, and Li fragments were measured 
and compared with both previously published measurements, and model predictions. This study 
is of importance for validating and improving carbon ion therapy treatment planning systems, and 
for estimating the radiological risks for personnel on space missions, since carbon is a significant 
component of the Galactic Cosmic Rays. 

PACS numbers: 25.70.Mn, 25.75.-q 

Keywords: Nuclear reaction, Nuclear emulsion, Heavy ion 
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I. INTRODUCTION 



To understand the fragmentation of heavy ions in high energetic nucleus-nucleus inter- 
actions is very important for many different applications, e.g., heavy ion radiotherapy 
and shielding against GCR (Galactic Cosmic Rays) ^nd SPE (Solar Particle Events) during 
manned missions outside the earth magnetic field [3']. For example, carbon ions of several 
hundred MeV/nucleon are currently used for radiotherapy at several facilities in the world. 
Radiotherapy with beams of heavy ions provides highly localized dose distributions at the 
end of the range (called Bragg peak), and is expected to be biologically more efficient than 
conventional photon radiation or proton beams. However, even if the heavy ions have the 
advantage of giving a localized dose distribution from the primary particles, the projectile 
fragmentation products can cause an undesired dose beyond the Bragg peak, which cause 
an unwanted dose to the healthy tissues. This effect can be estimated if the charge-changing 
cross sections of carbon ions in tissue equivalent materials, e.g., water, is known. Knowledge 
of fragmentation of carbon ions is also important when estimating the radiological risks in 
space, since the carbon is one of the most abundant heavy ion components in GCR, and 
its fragmentation products can make a significant contribution to the dose given to the 
personnel in the space vehicles. However, the fragmentation of carbon ion is not yet fully 
understood. Several measurements of the charge- changing cross sections of carbon in water 
and tissue materials have been performed 4, l5| . If compared to the predictions of different 



models 



0, y, y, 0, 



significant discrepancies are observed. Moreover, measurements 



of partial projectile fragmentation cross sections are scarce. More experimental data are 
therefore needed to better understand the fragmentation of carbon ion. 

From the experimental point of view, emulsion is one of the most suitable devices for de- 
tailed characterization of the fragmentation reactions, since it has the capability of spatial 
resolution of micron-order and excellent multi-particle separation in nearly Arc solid angle. 
Despite these advantages, this technology has not been applied in the past to a high statistic 
experiment because the track recognition is quite time-consuming due to the manual scan- 
ning of emulsion. However, the recent success of developing an automatic image processing 
system called Ultra Track Selector (UTS) 12] and an event reconstruction software package 
(NETSCAN) jl3| improved the scanning procedure greatly. Emulsion detectors are now used 
as core devices for large-scale elementary particle experiments 



fisj l. Taking advan- 
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tage of this new emulsion technique, we conducted an experiment to measure fragmentation 
of carbon ions at HIMAC (Heavy Ion Medical Accelerator in Chiba) at NIRS (National 
Institute of Radiological Sciences), Chiba, Japan. In this experiment, we constructed an 
emulsion detector, called Emulsion Cloud Chamber (ECC), which has the ability of three 
dimensional track and vertex reconstruction and charge identification for each individual 
track. 

In this paper, we report on measurements of total and partial charge-changing cross 
sections for incident ^^C beams in water and polycarbonate at energies ranging from 200 
to 400 MeV/nucleon. The results are discussed and compared with previous measurements 
and model predictions. 



II. EXPERIMENT 



A. Emulsion Cloud Chamber 



The Emulsion Cloud Chamber (ECC) was designed so that it can store information 
about all fragmentation reactions over the whole region where primary carbons penetrate 
the chamber. Schematic drawings of the ECC are shown in Fig. [H The ECC has 65 layers^ 
of emulsion modules and acrylic frames, and its total length is 29.1 cm long. The acrylic 
frame has a thickness of 2 mm and its inside was filled with de-ionized water at the exposure 
time. The ECC was installed in a water tank which has a window of 1-mm-thick acrylic 
plate on the side of beam injection. 

Each emulsion module has four emulsion sheets reinforced with a 1.05-mm-thick polycar- 
bonate (C16H14O3) plate whose measured density is 1.16 g/cw?. A pair of emulsion sheets 
covers on both sides of the plate as shown in Fig. [U^b).^ The polycarbonate plate keeps the 
flatness of the emulsion module, though it was also used as a target material. Each emulsion 
module is vacuum-packed with an aluminum-coated film for light and water shielding. The 



Iq . was 



OPERA film |l^ - a nuclear emulsion film developed for the OPERA experiment 
used as the emulsion sheet. An OPERA film measures 102 mm by 127 mm and has 44-/im- 
thick emulsion layers coated on both sides of a 205-/im-thick TAC (cellulose triacetate) base 



^ Consecutive numbers are assigned to the layers along the beam direction. 

^ The labels "A" to "D" are assigned to emulsion sheets in each module as shown in Fig.[ljb). 
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as shown in Fig. [U^c). 

B. Beam exposure and development 

The beam exposure was performed in the SB2 beam hne of the HIMAC heavy 
ion synchrotron in December 2004. The primary ^^C beam with a kinetic energy of 
400 MeV/nucleon was extracted to the beam hne and directly dehvered to the experimental 
area. 

The beam size and intensity were optimized so that the surface of ECC corresponding to 
the scanning area (12.5 cm^) was irradiated uniformly with a proper density. Beam density 
is one of the most critical parameters for the exposure, and exposures at excessive densities 
cause overlapping of particles, which makes the track reconstruction difficult. We controlled 
the beam density to be of the order of 10^ particles/cm^ which is the optimal density to 
observe fragmentation reactions at energies above 200 MeV/nucleon. We placed a 1-mm- 
thick plastic scintillator between the end of the beam pipe and the ECC, to monitor the 
number of injected charged particles. The total count of beam tracks was 23874, and the 
beam density was measured to be 1.6x10'^ particles/cm^. The purity of ^^C of the beam at 
the front of the ECC was estimated to be 99.5% by taking into account the interactions of 
^^C in the materials between the vacuum beam pipe and the ECC. 

After the exposure, the emulsion sheets were shipped to the scanning facility at Nagoya 
University. Sheets-A and C were developed normally and sheets-B and D were processed 
applying so-called the refreshing method [l^ before development. Refreshing is a technol- 
ogy to erase minimum ionizing tracks in emulsion by means of forced fading under a high 
temperature and high humidity environment. It desensitizes emulsion and extends the dy- 
namic range of response to highly charged particles. In Ref. jl8j] we showed that charges of 
ion tracks can be identified up to Z = 6 with the refreshing method. Sheets-B and D were 
refreshed for three days under the temperatures of 40 °C and 38 °C respectively and the 
relative humidity of 98%. 
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III. ANALYSIS 



A. Emulsion read-out 

We used the UTS (Ultra Track Selector) [12I] to read out track images in the emulsion 
films. UTS takes 16 layers of tomographic CCD images from a 44-/im-thick emulsion layer 
of a film. The imaging processor of UTS automatically finds track segments in an emulsion 
layer as three dimensional vectors. The UTS system has the field view of 150 fim by 120 /im 
and processes three views per second. Track segments with the tangents of the polar angles 
less than 0.5 were recorded, and the detection efficiency of a track segment was 98%. After 
the scanning, two corresponding track segments in the emulsion layers of both sides were 
connected across the TAC base, a new track segment was thus created. This is called a base 
track. The positional and angular resolutions of a base track were found to be 1 fim and 
6 mrad, respectively. 

UTS also records the grain density of a track segment, which keeps information on the 
local energy deposited. The grain density is digitized to pixel counts of CCD images whose 
pixel size is 0.29 /im by 0.23 /im. The pulse height is defined as the total pixel counts along 
each base track and used to identify a track charge at the later stage. 



B. Track and vertex reconstruction 

The NETSCAN software framework [l^ was utilized for track and vertex reconstruction. 
NETSCAN reconstructs trajectories of charged particles by connecting base tracks which 
lay within a certain distance. At the tracking stage, only base tracks in sheets-A were used. 
The probability of misconnection of two base tracks in adjacent modules was expected 
to be less than 1% for the primary beam tracks. The base tracks in sheets-B and D, 
which were desensitized with the refreshing method, were matched and associated with each 
reconstructed track. 

After the track reconstruction, the three dimensional vertex reconstruction was per- 
formed. Tracks starting inside the ECC and having a certain distance from a beam particle 
were selected as secondaries. Then, the vertex position was calculated searching for the 
position at which multiple tracks converged. We also imposed that an impact parameter 
of each secondary track should be less than 20 fim, which was determined to associate the 
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secondary tracks having the energy above 200 MeV/nucleon with a correct vertex at more 
than 99% efficiency by considering the effect of multiple coulomb scattering. The position 
resolution of a vertex was found to be 50 /um in the beam direction. 

C. Charge identification 

Pulse height information was used to identify the charge of each track. Figure [2] shows the 
pulse height distribution for secondary tracks, where the pulse heights are averaged among 
the connected base tracks in non-refreshing films of Sheets-A. As seen in the figure, pulse 
heights for multiple charged tracks overlap with each other and the track charges cannot be 
distinguished by this data alone. We therefore determined the track charges in two steps. 
At a ffist step, single charged tracks were discriminated from multiple charged tracks, de- 
manding that the averaged pulse heights are less than 210 as shown in Fig. [21 We estimated 
that 6% of single charged tracks were misidentified as multiple charged tracks by this sep- 
aration. Then, as a second step, we used pulse height information of refreshed films. The 
refreshing method desensitizes the emulsion films and decreases the pulse heights depending 
on refresh conditions. Sheets-B and sheets-D were refreshed in different conditions, so that 
they would show different responses to ion charges. Figure [3] presents a scatter plot between 
pulse heights averaged in sheets-B and sheets-D for secondary tracks. Clear discrimination 
of track charges from Z = 2 to Z = 6 can be obtained as seen in the figure. Likelihood 
analysis was applied for actual charge determination. The likelihood function is defined as 

L{Z)= n P^iZ,x), 

i=B,D 

where Pi{Z, x) is the probability function of a pulse height {x) for a specific charge [Z) and a 
refresh condition {i). We assumed that the pulse height distributions were Gaussian, where 
the mean pulse height and its deviation for individual charge were estimated by measured 
values for each refresh condition. Finally, the track charges were determined by varying the 
Z from two to six and finding the value which gave the maximum of the above likelihood 
function. The separation between Z = 5 and Z = 6 is most crucial, and the misidentification 
probability was estimated from the breadth of the pulse height distributions and found to 
be 1.4%. 
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D. Counting of charge-changing interactions 

The charge-changing interactions were identified principally by finding vertices whose 
primary tracks were observed as Z = 6 and any of its secondary tracks were not Z = 6. To 
ensure beam tracks, primary tracks were also required to be connected from the uppermost 
module. However a certain number of charge- changing interactions were missed using this 
identification condition. For instance, if a primary carbon and a secondary track were 
reconstructed as a single straight track, then the vertex would not be found. The following 
case is conceivable; a primary carbon changes to a secondary fragment with a kink angle 
too small to detect and no other secondary tracks are found in the angular acceptance 
(tan^ < 0.5). In order to recover these kinds of events, pulse height variations were traced 
along primary carbon tracks. Figure H] shows an example of this type of track (we call it 
a "charge-changing track"). As seen in the figure, a charge-changing interaction can be 
detected as a discrete change of pulse heights like a step function. It turned out that 87% 
of the interactions were detected as reconstructed vertices and the rest were counted as 
charge- changing tracks. 

For reconstructed vertices, we can recognize carbon-water and carbon-polycarbonate in- 
teractions by the vertex positions. Figure [5] shows the distribution of local distance between 
the surface of each module and reconstructed vertices (also see Fig. [H^b)). The target ma- 
terials of interactions were identified according to the detector structure as shown in the 
figure. For charge-changing tracks, we cannot directly determine the target material, be- 
cause we cannot get the exact vertex position. To evaluate the number of interactions, we 
did the following. The number of interactions, Nxii), in a target material of T (water or 
polycarbonate) of i-th module can be obtained by 

where N^'^ denotes the number of reconstructed vertices in a target (T) and N'^'^^ is the 
number of charge-changing tracks. Here Rt is introduced as the ratio of reactions in water 
or polycarbonate to those in the whole region of the ECC.'^ The ratio was given for each 

In Fig. [5l one can see a deterioration in vertex reconstruction at the first emulsion layer. This is due to 
interactions inside the emulsion layers which were used for tracking. Actual counts of interactions in the 
first emulsion region were taken as those in the second region. 
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charge- changing; = 1, 2, 3, 4, and 5. The counting inefficiency due to problems in 
reconstruction of secondary particles was estimated to be less than 2%. 

In the lower energy region, the track and vertex reconstruction deteriorates due to larger 
scattering. As effective events, we analyzed interactions in the upstream modules from 
35th module, which corresponds to the beam energy above 200 MeV/nucleon. Finally, we 
obtained 8213 interactions in this region of the ECC. 



IV. RESULTS AND DISCUSSION 

A. Total charge-changing cross section 

The total charge-changing cross section (atot) for a target material T in i-th module can 
be approximated by 

""'^'^'^ - N^pt NcitY 

where A^i is Avogadro's number, M is the molecular mass in atomic units of the target 
material, p is the target mass density, t is the target thickness, and Nc{i) is the number 
of primary carbon tracks impinging on i-th module. To reduce statistical errors, successive 
five (ten) modules are combined for carbon-water (carbon-polycarbonate) interactions. In 
Figs. E] and [7] our measured total charge-changing cross sections, from the interactions of 
carbon with water and polycarbonate, are plotted together with measurements from Ref. ^ 
and Ref. jsj. These cross sections and their statistical errors are also listed in Tables HI and HTl 
The cross sections of the carbon-polycarbonate (C-C16H14O3) reactions in the references are 
obtained from the combination of cross sections in different target materials; carbon-carbon 
(C-C), carbon-paraffin (C-CII2), and carbon- water (C-II2O) interactions. The beam energy 



at each module was estimated by energy loss calculation using Geant4 191]. As can be seen, 
our measurements agree well with the previous measurements at overlapping energies around 
250 MeV/nucleon. 

We have also compared our results with calculations based on the Sihver's model Q], 
which are superimposed in Figs. [6] and [71 Sihver's model is a semi-empirical model which 



takes advantage of the experimentally verified weak factorization property 2Q|]. The model 
is used in the treatment planning systems for heavy ion radiotherapy at several facilities in 
the world. It is also used in models for calculations of Galactic Cosmic Rays, etc. As can 
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be seen in Figs. |6] and [TJ the total charge-changing cross sections are quite consistent with 
the model predictions at energies above 250 MeV/nucleon, while the predictions differ from 
our results by about 10% in the lower energy region. 

The systematic error in our measurements consists of the following components: 

1. The uncertainty in vertex positions affected the position cut for selecting a target 
material, and resulted in an uncertainty of 3%. 

2. The probability of charge misidentification of primary or secondary particles gave the 
uncertainty of 5%. 

3. The inefficiency of counting charge-changing reactions was taken as the fractional 
systematic error of 2%. 

The contamination of and ^^C in the ^^C beam was expected to be 3%. Assuming 
that the differences of the cross sections between these isotopes are small, the contribution 



to the systematic error is considered to be negligibly small 
uncertainties, we get an estimated total systematic error of 6%. 



q]. By summing up all these 



B. Partial charge-changing cross section 

The partial projectile charge- changing cross section {(Taz) is defined as the cross section 
producing a projectile-like fragment with a specific charge difference AZ = 6 - ^fragment, 
ranging from AZ = 1 for B to AZ = 5 for protons. Figures [51fTU] and [TTlfT51 show our 
measurement results for carbon reactions of AZ = 1, 2, and 3 in targets of water and poly- 
carbonate, respectively. In Figs. [H [9], [11], and [121 our measured partial charge-changing 
CO. .ct.o,. a. *o co,npa.ed w.th the .e.u.e.e„t. pubLshed . Ref. g AH ™. 
sections are also listed in Tables [1111 and [IVl As can be seen from these figures, our measure- 
ments agree well with previously performed measurements at overlapping energies around 
250 MeV/nucleon. 



We have also compared with calculations based on the Sihver's model 



6], which are 



superimposed in Figs. [8][T3l^ The model predictions for the partial cross sections show 



^ The values of Sihver's model were obtained by summing up for all detectable isotopes with A = 11, 10, 8 
for AZ = 1. A = 10, 9, 7 for AZ ^ 2. and A = 7. 6 for AZ = 3. 
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stronger energy dependence than our measurements, and we found significant discrepancies, 
especially in the lower energy region. The deviations between the measured and the calcu- 
lated cross sections for the fragments with AZ = 1, 2, and 3 at the projectile energy 255 
MeV/nucleon, are 21%, 23%, and 24% for the carbon- water interactions, and 16%, 33%, 
and 4% for carbon-polycarbonate interactions, respectively. The total systematic error was 
estimated to be 6%, in the same way for the total cross sections. 

We do not present the results of the partial charge- changing cross sections for AZ = 4 
and 5 i.e., the production of a and p, since these light particles have energy and angular 
distributions beyond the detector acceptance; and contamination of light target fragments 
can not be neglected. 

V. CONCLUSIONS 

We have measured charged nuclear fragments produced by 200 - 400 MeV/nucleon carbon 
ions, interacting with water and polycarbonate, using a newly developed emulsion technique. 
Both total and partial charge-changing cross sections were measured, and the results were 
compared with previous measurements and model predictions. The measured total charge- 
changing cross sections showed good agreement with both the reference measurements and 
simulations using a model developed by Sihver et al. [6]. In the case of the partial cross 
sections Sihver's model showed stronger energy dependence than our measurements, and 
deviated from our measurements in the lower energy region. This proves that measurements 
performed by our ECC can help improve and validate theoretical fragmentation models. 
Since 200 - 400 MeV/nucleon is the energy region where carbon ion therapy is performed 
at several facilities, our data is of importance for validating and improving the treatment 
planning systems. The results are also of importance when estimating the radiological risks 
for personnel on space missions, since carbon is a significant component of the GCR. 
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TABLE I: Comparison of total charge-changing cross sections of C beams in water target. 



Our measurements 


Ref. [4] 




Ref. [5j 




E (MeV/nucleon) 


a tot (mb) 


E (MeV/nucleon) 


dtot (mb) 


E (MeV/nucleon) 


atot (mb) 


377 ± 12 


1253 ± 45 










352 ± 13 


1202 ± 47 










326 ± 13 


1250 ± 51 










299 ± 14 


1183 ± 52 










271 ± 15 


1193 ± 56 


259 ± 47 


1205 ± 34 






241 ± 16 


1241 ± 61 


241 ± 26 


1163 ± 13 


241 ± 33 


1191 ± 35 


208 ± 17 


1231 ± 64 











TABLE IL Comparison of total charge-changing cross sections of ^^C beams in polycarbonate 
(C16H14O3) target. The values of Ref. [4 R-^f. 5] are obtained from the combination of cross 
sections in different target materials (C, CH2, and H2O). 



Our measurements 


Ref. [4] 


Ref. [5J 


E (MeV/nucleon) 


(Jtot (mb) 


E (MeV/nucleon) atot (mb) 


E (MeV/nucleon) atot (mb) 


364 ± 25 


17682 ± 595 






312 ± 27 


16723 ± 648 






255 ± 30 


17225 ± 747 


226 ± 41 16765 ± 134 


241 ± 33 17625 ± 562 
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TABLE III: Comparison of partial charge-changing cross sections of C beams in water target. 



Om' measurements Ref. [5] 



AZ E (MeV/nucleon) 


CTAz (mb) 


E (MeV/nucleon) cjaz (mb) 


1 364 ± 25 


221 ± 14 




312 ± 27 


218 ± 15 




255 ± 30 


210 ± 17 


241 ± 33 210 ± 6 


2 364 ± 25 


89 ± 9 




312 ± 27 


95 ± 10 




255 ± 30 


93 ± 11 


241 ± 33 84 ± 4 


3 364 ± 25 


120 ± 10 




312 ± 27 


115 ± 11 




255 ± 30 


128 ± 13 





TABLE IV: Comparison of partial charge-changing cross sections of ^■^C beams in polycarbonate 



(C16II14O3) target. The values of Ref. 



5(] are obtained from the combination of cross sections in 



different target materials (C, CII2, and II2O). 



Our measurements Ref. [5] 



AZ E (MeV/nucleon) 


fJAZ (mb) 


E (MeV/nucleon) 


<TAZ (mb) 


1 364 ± 25 


2968 ± 243 






312 ± 27 


3162 ± 281 






255 ± 30 


2714 ± 296 


241 ± 33 


2890 ± 70 


2 364 ± 25 


957 ± 138 






312 ± 27 


1158 ± 170 






255 ± 30 


1027 ± 182 


241 ± 33 


1028 ± 43 


3 364 ± 25 


1903 ± 195 






312 ± 27 


1798 ± 212 






255 ± 30 


1992 ± 253 
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(a) ECC whole structure 
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(b) Detailed structure 
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FIG. 1: (Color online) (a) A schematic view of the whole structure of the Emulsion Cloud Chamber 
(ECC). It has 65 layers of emulsion modules and water targets, (b) Detailed structure of the ECC. 
An emulsion module consists of four emulsion sheets and a polycarbonate plate, (c) Cross section 
of an OPERA film. An OPERA film has emulsion layers coated on both sides of a TAC base. 
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FIG. 2: Distribution of the pulse heights in sheets- A for secondary tracks. The tracks with pulse 
heights less than 210 are discriminated as single charged tracks. 
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FIG. 3: A scatter plot between averaged pulse heights in sheets-B and sheets-D for secondary 
tracks. One can see clear discrimination of track charges from Z = 2 to Z = 6. 
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FIG. 4: An example of a charge-changing track. Open circles (Black squares) represent averaged 
pulse heights in sheets-B (sheets-D). We traced pulse height variations along primary carbon tracks. 
In this carbon beam changes into a beryllium between 24th and 25th modules. 
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FIG. 5: Distribution of local distance between the surface of each module and reconstructed 
vertices. The target materials of interactions were identified according to the detector structure 
(also see Fig. [U^b)). 
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FIG. 6: Total charge-changing cross sections of carbon-water interactions in comparison with 

previous measurements and a model calculation. The vertical error bars indicate statistical errors 
only in this and the followings figures. 
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FIG. 7: Total charge-changing cross sections of carbon-polycarbonate interactions in comparison 
with previous measurements and a model calculation. 
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FIG. 8: Partial charge-changing cross sections of carbon-water interactions of AZ = 1 in comparison 
with previous measurements and a model calculation. 
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FIG. 9: Partial charge-changing cross sections of carbon-water interactions of = 2 in comparison 
with previous measurements and a model calculation. 
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FIG. 10: Partial charge-changing cross sections of carbon- water interactions of A = 3 in com- 
parison with a model calculation. 
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FIG. 11: Partial charge-changing cross sections of carbon-polycarbonate interactions of AZ = 1 
in comparison with previous measurements and a model calculation. 
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FIG. 12: Partial charge-changing cross sections of carbon-polycarbonate interactions of AZ = 2 
in comparison with previous measurements and a model calculation. 
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FIG. 13: Partial charge-changing cross sections of carbon-polycarbonate interactions of AZ = 3 
in comparison with a model calculation. 
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